Large-scale anisotropy in scalar turbulence 
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The effect of anisotropy on the statistics of a passive tracer transported by a turbulent flow 
is investigated. We show that under broad conditions an arbitrarily small amount of anisotropy 
propagates to the large scales where it eventually dominates the structure of the concentration field. 
This result is obtained analytically in the framework of an exactly solvable model and confirmed by 
numerical simulations of scalar transport in two-dimensional turbulence. 

PACS numbers: 47.27.-i 



The emergence of large-scale anisotropy arising from 
small-scale sources is a phenomenon that spans the most 
diverse fields of physics. For instance, the microscopic 
anisotropy of crystals in mantle rocks in the Earth's in- 
terior is believed to induce large-scale seismic anisotropy 
0, and some small anisotropic perturbation in an early 
cosmological era, evolving through gravitational collapse, 
is thought to be responsible for the large-scale anisotropy 
of the cosmic microwave background radiation 0] ■ Con- 
versely, in statistical physics, microscopic details such as 
lattice anisotropics may be wiped out by the dynamics al- 
lowing a recovery of symmetry and universality at large 
scales In the theoretical and experimental analysis 
of turbulence much attention has been devoted to the 
anisotropy of the fine scales of fluid motion (see, e.g., Q 
and references therein). Here, we take a different view- 
point and investigate the effect of anisotropy on the large- 
scale statistics of turbulence. In this Letter we show how, 
unexpectedly, breaking rotational invariance by an arbi- 
trarily small amount at a given scale induces a strong 
anisotropy on the large scales, and symmetry is never 
restored. 

We consider the evolution of a passive tracer described 
by a concentration field 6(x,t) and transported by a tur- 
bulent flow v(x,t) 



d t 9 + v- Ve = KA0 + f, 



(1) 



where v is an incompressible, statistically homogeneous 
and isotropic velocity field. The external driving / is the 
source of scalar field fluctuations acting at a character- 
istic scale I f . The turbulent cascade toward small scales 
produces fine-scale structures of concentration that are 
eventually smeared out by diffusion at scales r& <C If, 
resulting in a statistically stationary state where input 
and dissipation are in balance on average. The pumping 
mechanism can be chosen so as to introduce a certain 
degree of anisotropy, that propagates across scales and 
may in principle pervade the system. However, the dis- 
ordered motion of fluid particles induced by the under- 
lying turbulent, isotropic medium, might be sufficient to 
restore rotational invariance at scales far below or above 



If. Indeed, this is the case at small scales r -C If, where 
it can be shown that the dominant contribution to the 
statistics of the scalar field 6 is isotropic At large 
scales r ^> If, since no upscale cascade of scalar fluctu- 
ations occurs, a fortiori one would expect an essentially 
isotropic concentration field. On the contrary, here we 
give theoretical and numerical evidence that large-scale 
statistics is dominated by the anisotropic contribution 
under very broad conditions. We show that: (i) The 
correlation function (0(r,t)9(0,t)) at scales r ^> If is 
dominated by its anisotropic component decaying as a 
power law with an anomalous scaling exponent, as op- 
posed to the exponential fall-off of the isotropic part. 
This result is obtained analytically in the framework of 
the exactly solvable Kraichnan model and its validity for 
realistic flows is demonstrated by numerical simulations 
of passive scalar advection in the inverse cascade of two- 
dimensional turbulence, (ii) Large-scale anisotropy man- 
ifests itself in the concentration field with the appearance 
of "pearl necklace" structures aligned with the prefer- 
ential direction imposed by the microscopic anisotropy. 
(Hi) The loss of isotropy at large scales can be interpreted 
as a breakdown of equilibrium Gibbs statistics for the 
anisotropic degrees of freedom; (iv) In the Lagrangian 
interpretation of passive scalar transport the emergence 
of anisotropic power-law decay of correlation is associ- 
ated to a long-lasting memory of the initial orientation 
of particle pairs advected by the flow. 

Let us first consider the Kraichnan model of pas- 
sive scalar advection (see, e.g., U for a thorough re- 
view), where v is a Gaussian, self-similar, incompress- 
ible, statistically homogeneous and isotropic, white-in- 
time, d— dimensional velocity field. Its statistics is char- 
acterized by the correlation S a p(r)5(t) = {[v a (r,t) — 
v a (O,t)][v (r,O) - w/j(0,0)]> = 2Dri[(d + £ - l)5 aP - 
^T a rp/r ]S(t). The exponent £ measures the degree 
of roughness of the velocity field and lies in the range 
< £ < 2, the two extremes corresponding to Brownian 
diffusion and smooth velocity, respectively. The assump- 
tion of ^-correlation is of course far from being realistic, 
yet it has the remarkable advantage of leading to closed 
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equations for equal-time correlation functions. In the 
following it will be sufficient to focus on the two-point 
correlation function C(r) = (9(r,t)6(0,t)). In the limit 
of vanishing diffusivity n — > and in the statistically sta- 
tionary state, C obeys the partial differential equation 
MC(r) = -F(r) where M = \S a(3 (r) Here F is 
the correlation function of the Gaussian, white-in-time, 
statistically stationary, homogeneous, anisotropic forcing 
(f(r,t)f (0,0)) = F(r)S(t). At scales r < l f it equals the 
average input rate of scalar and then decays rapidly to 
zero, e.g. exponentially, as r > lj. By virtue of the 
statistical isotropy of the velocity field, the operator M 
assumes a particularly simple form in radial coordinates: 
M = D[(d-l)r 1 - d d r r d -' l - + Sd r + (d+£-l)rt- 2 L 2 ] where 
L 2 is the d-dimensional squared angular momentum op- 
erator. It is then convenient to decompose the corre- 
lation functions on a basis of eigenfunctions of angular 
momentum £-> 2 Yj — —j(j + d — 2)Yj with positive integer 
j. The short-hand notation Yj(r), where r = r/r, does 
not account for degeneracies and stands for the trigono- 
metric functions in d = 2 and the spherical harmonics 
in d — 3. Accordingly, we define the components of 
the correlation functions in the j-th anisotropic sector 
as C(r) = ^ZjCj(r)Yj(r) and similarly for F(r), where 
Cj and Fj depend on r = \r \ only. This yields a system of 
uncoupled differential equations in the radial variable for 
each anisotropic component JAjCj(r) = —Fj(r), where 
= D[(d - iy-dd.r4-i+£^. _ j(j + d _ 2 )(d + £ - 
l)r^~ 2 ], that can be solved in each sector j. The result- 
ing Cj (r) is a linear combination of a particular solution 
determined by Fj(r), and a homogeneous one Zj(r), a 
"zero mode". It is easy to see that the former behaves 
as r 2 ^^ +J for r < (/ (recall that Fj ~ r 3 at small r if 
F is analytic in the neighborhood of r = 0) and that 
it must fall off exponentially fast for r ^> If , as dic- 
tated by the decay of Fj. The homogeneous solutions 
are Z^(r) = with scaling exponents Cf — 5!^ + 
M+d-Wd+t-i) , The Z ero mode 



2-£±^(d-2 + £) 2 +4 
with positive scaling exponent appears at small scales 
whereas the zero mode with negative scaling exponent C~ 
is relevant in the range r^>lf. In order to fully char- 
acterize the large-scale behavior of the correlation func- 
tion, it is necessary to identify the prefactor appearing in 
front of the homogeneous solution. This can be accom- 
plished by writing the equation for Cj in integral form: 
Cj(r) = J °° Gj(r, p)Fj(p)dp, where Gj(r,p) is the kernel 
of —MJ 1 , i.e. the solution of M,jGj(r, p) = —5(r — p). 
The explicit form is Gj(r, p) — A(p) Z J + (r) Z J _ (p) for r < p 
and Gj(r,p) = A(p)Zi{r)Z 3 + {p) for r > p, with A(p) = 
p d ~ l /[£ ) (C J + — Cj 7 )]- Plugging this expression in the in- 
tegral form for the correlation function yields a large- 
scale behavior Cj(r) ps QjZJ(r) + exponentially decay- 
ing terms. The quantity Qj — J °° A(p)Z 3 + (p)Fj(p) dp is 
of crucial importance and plays the role of a "charge" 




FIG. 1: Image of a scalar field for the Kraichnan model at 
£ = 0, d — 2, Qo = 0, Q2 7^ 0. The width of a single scalar 
patch is ~ If. Forcing is preferentially acting in the vertical 

direction. Here Co(r) ~ e~ r and C%(r) ~ r~ 2 at r > I/. 



in analogy with electrostatics |6| . In the isotropic sector 



r(rf/2) 



■ / F(r) dr. When 



j = 0, it reduces to Q = 2wd/ 2 D(d _ 2+a 
the isotropic charge Qo ^ 0, the leading behavior at large 
scales is isotropic, C(r) ~ Co(r) ~ Qor~ d+2 ~^ . The 
most interesting situation is when Qq — 0, correspond- 
ing to the broad class of forcings localized in wavenum- 
ber space (Qo °c F(k = 0)). In this event, there is no 
power-law contribution from the isotropic zero mode and 
therefore the isotropic part of the correlation function is 
characterized by an exponential decay at large r 0,0- In 
the anisotropic sectors, it appears immediately that there 
is no reason to expect a null charge and the generic situ- 
ation is Qj 7^ for j > (see 8]). As a result, the large- 
scale correlation is dominated by the anisotropic contri- 
bution arising from the zero- mode, C(r) ~ r^j Yj(r), 
that largely overweights the exponentially small isotropic 
part. Among the various contributions arising from dif- 
ferent sectors, the leading one corresponds to the lowest 
nonzero j excited by the forcing, typically j — 2 (odd j's 
are switched on only by breaking reflection invariance). 

For the sake of illustration, we show in Fig. ^ an in- 
stance of a scalar field corresponding to the simple case 
of Kraichnan advection by a very rough velocity (£ = 0), 
in d = 2 and with Qo = 0. Large-scale anisotropy man- 
ifests itself in the appearance of "pearl necklaces" made 
of like-sign scalar patches of size ~ If. These are aligned 
along the preferred direction of the forcing and extend 
for a length ^> If . 

It is worth pointing out the relationship between 
the appearance of anisotropic, anomalous scaling in the 
large-scale behavior of the scalar correlation function and 
equilibrium statistics. At large scales, because of the ab- 
sence of scalar flux, the system could be expected to be in 



3 



equilibrium and obey Gibbs statistics. In physical space 
this corresponds to a concentration field organized in in- 
dependently distributed scalar patches of size If. As 
recently shown in Refs. 0, Q , this is not true and sub- 
stantial deviations are observed at the level of multipoint 
correlation functions already in the isotropic case. This 
departure has been traced back to the existence of non- 
trivial zero-modes in that case as well. In the present 
case the breakdown of Gibbs statistics has an even more 
dramatic manifestation since it occurs already for two- 
point correlation functions, i.e. at the level of the spec- 
tral distribution of concentration. We now rephrase the 
previous findings in terms of the averaged scalar spec- 
tral density, i.e. the Fourier transform of the correla- 
tion function, C(k) — (\0(k, i)| 2 ), and its decomposition 
in angular sectors in wavenumber space ■ Cj(fc)Y} (fc), 

where k = k/k. For a correlation function decaying ex- 
ponentially to zero at large r, representative of large- 
scale equipartition in physical space, the spectral density 
is analytic in a neighborhood of k = 0. In the series 
for C(k) the harmonic Yj(k) appears only in the powers 
of k of order > j, yielding the long wavelength behav- 
ior Cj eq \k) ~ k J . This defines the equipartition spec- 
trum for generic anisotropic fluctuations. However, be- 
cause of the appearance of nontrivial zero modes in the 
anisotropic sectors the actual spectral density contains 
also a contribution Cj Ze (k) ~ fc - ^^ that is respon- 
sible for the power-law decay of correlations in physical 
space with j > 0. For the Kolmogorov-Richardson value 
£ = 4/3 and the sector j = 2 the anomalous spectrum 
always dominates the equipartition contribution in spec- 
tral space as well (— d — CJ=2 < 2 for all £ < 3/2). 

It is useful to reinterpret the results obtained so far 
within the framework of the Lagrangian interpretation 
of passive scalar transport. The correlation function can 
be generically written as C(r) = J T(p\r) F(p) dp where 
T(p\r) dp is the average time spent at a separation p+dp 
by a pair of particles that end their trajectories at a 
separation r. In the Kraichnan model T is the kernel 
of the operator — M" 1 . Since the action of the forcing 
is restricted to scales ~ If, the large-scale behavior of 
the correlation function is essentially dominated by the 
ensemble of trajectories that have spent in the past a 
sufficiently long time at a short distance p < If <C r. 
When r/p tends to infinity T becomes independent of 
p and we obtain T ~ r ~ d+2 ~^ The dependency on 
the final orientation r is also lost. This leads to the 
estimate C(r) ~ Qor~ d+2 ~^ . However, as noticed previ- 
ously, when Qq = the isotropic part of the correlation 
function receives only exponentially small contributions 
from the forcing in the range p^r^>lf. Let us now 
turn our attention to the anisotropic part of the corre- 
lation function. Projecting C(r) over Yj(f) for j > 
amounts to give different weigths, positive and negative, 
to particle pairs oriented in different directions r. Thcrc- 




10~ 2 10" 1 10° 

kl, 

FIG. 2: Simulations of passive scalar advection by two- 
dimensional Navier-Stokes turbulence in the inverse cascade 
(for details see 0). Here are shown the spectra Ej(k) of con- 
centration in the isotropic (j — 0) and anisotropic (j = 2) 
sectors. The spectral slopes are 1.0 ± 0.05 and 2.2 ± 0.1, re- 
spectively (Cj^2 = — 3.2). The equilibrium spectra k and fc 3 
are also shown for comparison. 



fore Cj can be interpreted as a difference of times spent 
at p < If by differently oriented pairs. The first key point 
is that the trajectories preserve a long-lasting memory of 
their initial orientation, with a slow power-law decay in 
r that reflects in the behavior of the correlation function. 
Indeed, it can be shown that in the Kraichnan model 

T{p\ r ) = p c i Yj(f)Yj(p) for r > p. Plugging 

this expression in the integral form of the correlation 
function gives the result Cj(r) ~ Qjr^i as above. Here 
emerges the second important point, i.e. the dependence 
of T on p: differently oriented trajectories sample the 
forcing unevenly in scales as p ~ If and this results in a 
nonvanishing charge Qj for j > 0. 

A remarkable advantage of the Lagrangian interpreta- 
tion is that it does not make appeal to the special features 
of the Kraichnan model. This suggests that the same 
mechanisms are at work for realistic turbulent flows as 
well, and this expectation has been repeatedly confirmed 
for different aspects of passive scalar transport 0,EJE1- 
Here we show that anisotropy dominates the large-scale 
statistics for real flows by showing the results of a nu- 
merical investigation of passive scalar transport in the 
inverse cascade of two-dimensional Navier-Stokes turbu- 
lence. This flow has been studied in great detail, both 
experimentally in fast flowing soap films ^3] and hi shal- 
low layers of electromagnetically driven electrolyte solu- 
tions Hilllil], and numerically 0EEII3- The velocity 
field v is statistically homogeneous, isotropic, and scale- 
invariant with exponent ft, = 1/3 (6 r v ~ r ) in the range 
~ r ^ where IJ denotes the kinetic energy in- 
jection length and L v the velocity integral scale. The 
scalar field is governed by Eq. Q driven by a homoge- 
neous, anisotropic, Gaussian, (^-correlated driving / that 
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In summary, we have shown that microscopic 
anisotropics introduced by the forcing have a dramatic 
imprint on the large-scale statistics of passive scalar tur- 
bulence. From this result new questions arise naturally, 
the most intriguing one being whether the large scales of 
hydrodynamic turbulence show such striking properties 
as well. Further theoretical, experimental and numerical 
effort in this direction is needed to elucidate this point. 

This work has been partially supported by the EU con- 
tract HPRN-CT-2002-00300 and by CNRS PICS 3057. 



FIG. 3: Main frame: Correlation functions Cj(r) at large 
scales r > If for j = and j = 2. Notice the exponen- 
tial decay of the isotropic part as opposed to the power-law 
behavior of the anisotropic component. In order to limit 
the statistical noise at large values of r/lf these two curves 
have been obtained by interpolating the spectra shown in 

Fig. |21 at small k as power laws k and fc -1 i= 2 , respec- 
tively, and computing the correlation function from the inte- 
gral Cj{r) = Jj{kr)Ej{k)/kdk. Here C~= 2 = "3.2. Inset: 
Correlation function Cj= 2 {r) at small scales r <lj. The solid 

line is r C J= 2 with £t 2 = 1.8. 



excites the sectors j = and j = 2 and satisfies the 
condition of null isotropic charge (see [§]). The various 
lengthscales are ordered as follows: <C Td <C // -C L v . 
In Fig. we show the spectral content of scalar fluctu- 
ations Ej(k) = n~ 1 k J Q * cos(j<fik) C(k) d<j)k = kCj(k) at 
klf < 1, i.e. at large scales. The isotropic spectrum 
(J = 0) agrees very well with the Gibbs equilibrium dis- 
tribution, Ej =0 (k) oc k, and corresponds to exponentially 
decreasing isotropic correlation at large scales (see the 
main frame of Fig. |3J) in agreement with the theoretical 
arguments presented above. The anisotropic spectrum 
shows a power-law behavior Ej = 2(k) ~ fc 2 - 2±01 definitely 
dominating over the equilibrium spectrum EjL 2 (k) ~ 
In physical space, this translates in a power-law decay of 
the anisotropic correlation function Cj—2(r), as shown 
in the main frame of Fig. and leads to the estimate 
Qj—2 — —3.2 ± 0.1. Therefore, the correlation function 
at large scales is dominated by the anisotropic power- 
law decay for 2D Navier-Stokes advection as well. Fi- 
nally, we notice that for incompressible, time-reversible, 
self-similar flows the two zero-mode exponents are conju- 
gated by the dimensional relation (j~ + C,J = — d + 1 — h 
(within the Kraichnan model 1 — h — 2 — £, due to the 
(^-correlation in time). In the inset of Fig. |21 we show 



the behavior of the correlation function C J= 2 (r) 



at small scales r < 1/ that indeed displays an exponent 
Cj~ = 1.8 ±0.1 compatible with the previous relation [l8| . 



[1] P. G. Silver, Ann. Rev. Earth. Planet. Set. 24, 385 (1996) 
[2] E. F. Bunn, P. G. Ferreira, J. Silk, Phys. Rev. Lett. 77, 
N 14 (1996) 

[3] J. Cardy, Scaling and renormalization in statistical 
physics, Cambridge University Press, Cambridge UK 
(1996). 

[4] L. Biferale and I. Procaccia, Phys. Rep. 414, 43 (2005) 
[5] G. Falkovich, K. Gawedzki and M. Vergassola Rev. Mod. 

Phys. 73, 913 (2001) 
[6] G. Falkovich and A. Fouxon, Phys. Rev. Lett. 94, 214502 

(2005) 

[7] A. Celani and A. Seminara, Phys. Rev. Lett. 94, 214503 
(2005) 

[8] A simple example is given by the following choice of 
anisotropic driving in d — 2: the spectral density of the 
forcing F(k,(j)k) = 7r _1 fc 2 ^ exp(— fc 2 Z^/2) sin 2 (f)k results 
in a forcing correlation function that excites the two 
sectors j = 0, 2 only: F(r,4>) = F (r) + F 2 (r) cos(2<?i), 
where F (r) = 2[1 - r 2 / (2l 2 f )] exp[-r 2 /(2l 2 f )} and 
F 2 {r) = -(r/; / ) 2 exp[-r 2 /(2^)]. The result- 
ing charges are Qq = (_D£) _1 L 00 fo{r)r dr oc 
J °°r[l - r 2 /(2^.)]exp[-r 2 /(2Z 2 .)]dr' = and 

Q 2 = (v / 2«/) 2+c ?r(2 + C 2 + /2)/[£(C 2 + - C a ")] + 0. 
[9] D. Bernard, K. Gaw§dzki, and A. Kupiainen, J. Stat. 

Phys. 90, 519 (1998) 
[10] A. Celani, A. Lanotte, A. Mazzino, and M. Vergassola, 

Phys. Rev. Lett. 84, 2385 (2000) 
[11] A. Celani and M. Vergassola, Phys. Rev. Lett. 86, 424 

(2001) 

[12] H. Kellay and W. Goldburg, Rep. Prog. Phys. 65, 845 
(2002) 

[13] P. Tabeling, Phys. Rep. 362, 1 (2002) 

[14] M. K. Rivera and R. E. Ecke, Phys. Rev. Lett. 95, 194503 

(2005); S. Chen et at, ibid. 96, 084502 (2006) 
[15] L. Smith and V. Yakhot, Phys. Rev. Lett. 71, 352 (1993) 
[16] G. Boffetta, A. Celani and M. Vergassola, Phys. Rev. E 

61, R29 (2000) 

[17] D. Bernard, G. Boffetta, A. Celani and G. Falkovich, 
Nature Physics 2, 124 (2006) 

[18] Notice that, strictly speaking, the assumption of time- 
reversibility for the inverse cascade is incorrect, due 
to the presence of a kinetic energy flux. However, the 
almost-Gaussian statistics of velocity (see, e.g., Ref. (l^l - ) 
suggests that it is a good approximation. 



